Robust to impurity-scattering spin Hall effect in two-dimensional electron gas 
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We propose a mechanism of spin Hall effect in two-dimensional electron gas with spatially random 
Rashba spin-orbit interaction. The calculations based on the Kubo formalism and kinetic equation 
show that in contrast to the constant spin-orbit coupling, spin Hall conductivity in the random 
spin-orbit field is not totally suppressed by the potential impurity scattering. Even if the regular 
contribution is removed by the vertex corrections, the terms we consider, remain. Therefore, the 
intrinsic spin-Hall effect exists being, however, non-universal. 
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Spin currents are believed to be of great importance 
for future spin electronics^, as they offer the possibil- 
ity of nonmagnetic manipulation of magnetic moments. 
Generally, spin currents are associated with charge cur- 
rents, and can be generated by various methods, such as 
for instance by electric field in magnetic systems or cir- 
cularly polarized light in nonmagnetic semiconductors. 
Of particular interest, however, are pure spin currents, 
where the flow of spins is not accompanied by any elec- 
tric current. Search for generation techniques of pure 
spin currents, especially in nonmagnetic semiconducting 
systems, is of high interest both for fundamental and ap- 
plied physics. One of the possibilities of producing pure 
spin currents relies on the spin Hall effect (SHE) in non- 
magnetic semiconductors with various types of spin-orbit 
(SO) interaction, where uniform electric field causes a 
transverse spin rather than charge current. 

The existence of SHE in semiconductors with impuri- 
ties has been predicted by Dyakonov and Perel— . Since 
then several mechanisms of SHE have been proposed^—, 
and the effect has been observed in a number of 
experiments^. It is generally believed that two extrinsic 
mechanisms related to the SO-dcpcndcnt scattering by 
impurities, i.e. side jump and skew scattering^—, can be 
responsible for the SHE in metallic and semiconducting 
materials. In addition, a lot of discussions in recent liter- 
ature concerned the possibility of SHE due to intrinsic SO 
interaction in disorder-free systems. An extensively stud- 
ied example of such a system is a two-dimensional elec- 
tron gas with constant Rashba SO interaction^ 1 ^, lead- 
ing to the momentum-dependent spin splitting of electron 
states. The theoretical efforts were especially focused on 
the possibility of equilibrium spin current s 13 ' 14 and uni- 
versal SHE£ independent on the SO coupling strength. 
However, it turned out that the role of impurities is cru- 
cial for this mechanism 1 ^. It has been shown that even in 
the limit of a very weak spin-independent disorder, the 
potential scattering from impurities suppresses the SHE 
completely-^— . In the case of random Rashba field with- 



out impurity scattering, the SHE can be also nonzero, 
as shown numerically within tight-binding model for a 
finite-size system^ 4 -. 

Here we show that the intrinsic SHE does exist, para- 
doxically, in relatively dirty systems, where the SO cou- 
pling appears locally, but vanishes on average. Such a 
random spin dynamics is common in symmetric semi- 
conductor quantum wells (QWs)2^, such as Si/SiGe2& 
and GaAs/AlGaAs QWs grown along the [110] axis^ 7 -. 
Moreover, we show that impurities play here the role less 
important than in the case of uniform SO interaction - 
the spin Hall conductivity does not vanish in the limit 
of small impurity density. This behavior is qualitatively 
different from that for constant Rashba SO interaction. 

To describe the model under consideration we assume 
Hamiltonian of electrons moving in the r = (x, y) plane 
with random Rashba SO interaction, H = Hq + H so (in 
the following we use units with h = 1), with 
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2m' 



(1) 



lcr x {V y ,X(v)}+ l -a y {V x ,X(r)}, (2) 



where Vj = <9; — eAj/c, A is the vector potential of ex- 
ternal field, e and m is the electron charge and effective 
mass, respectively, and a a are the Pauli matrices (a = 
x,y,z). The curly brackets {...} stand for the anticom- 
mutator of the appropriate operators to ensure the Her- 
mitian form of the Hamiltonian. The random coupling 
parameter A(r) vanishes on average, (A(r)) = 0, while 
the correlator C A a (r-r') = (A(r)A(r')) = (A 2 ) F(r - r'), 
with all higher correlators reduced to the second-order 
one for the Gaussian fluctuations of A(r). 

The spin current operator has the following form 2 ^: 



Ji = J~ {Hi, °a}, 



(3) 



where ji = —c(dH/dAi) is the i-th component of the 
current operator (i — x,y,z). We consider in-plane 
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electric field E, and calculate the total spin current J". 
In the following we use the gauge with vector potential 
A(t) = A e- iwt , E = -c -1 (9A/0i), and at the end take 
the limit uj —¥ in the calculated response function^. 

Using Eqs. ([T]) and ([3]) one can write the matrix ele- 
ments of the spin current operator in the basis of eigen- 
functions of Hq as 



(k|j*|k') = |f 



1 c 



(4) 



where k includes the electron momentum k and spin com- 
ponent <j z . We note that the z-component of spin cur- 
rent. J?, does not contain any anomalous part explicitly 
dependent on the SO coupling. 

It is convenient to decompose the Hamiltonian H into 
two terms, H = Ha=o+Ha, where Ha=o corresponds to 
vanishing vector potential (A = 0), while Ha appears at 
nonzero A. Matrix elements of Hamiltonian Ha=o are 
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k|J?A=o|k 

Akk 



k 2 



(5) 



kk' 



- [a x (k y + k' y ) - a v (k x + k' x )] , (6) 



where Akk' is the Fourier component of the ran- 
dom Rashba field. In turn, matrix elements of the 
A— dependent term, Ha, have the form 



A 



mc 
eAkk' 



u kk' ^ 

(0 X Ay 



2mc 2 kk 

~ &yA x ) . 



(7) 



(8) 



In the linear response regime, only the first and third 
terms on the right-hand side of Eq. ([7]) are relevant. The 
third term clearly demonstrates the coupling of electric 
field to electron spin via the Fourier component of the 
random Rashba field. These two terms can be associated 
with two different electromagnetic vertices in the Feyn- 
man diagrams for system's conductivity: the first one 
leads to the conventional conductivity while the third 
one to the spin conductivity. 

To calculate the spin Hall conductivity we apply the 
conventional Kubo formalism 2 ^ using the retarded and 
advanced Green's functions G 
vicinity of the Fermi level, 



R.A 



IG^' A , taken in the 



qR,A _ 1 

k e k -eF±i/2r' 



(9) 



where / is the 2x2 unit matrix, e k = k 2 /2m, ef is 
the Fermi energy, and r is the total momentum relax- 
ation time including scattering from impurities (to) and 
scattering by spin-dependent Rashba potential (rso)j 
1/t = 1/tq + 1/tso- Since the SO coupling vanishes 
on the average, the Green function © keeps exactly the 
diagonal form in the spin subspace. The linear spin Hall 
conductivity is represented by the sum of two Feynman 
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FIG. 1: The Feynman diagrams leading to nonvanishing con- 
tributions to spin current. Here the left vertex (filled square) 
corresponds to the spin current operator, the right vertex 
(filled circle) is the external field perturbation in Eq.©, and 
the white circle is the matrix element of spin-orbit coupling 
m Eq.©. Upon averaging over disorder, the dashed line 
becomes the Fourier component of the correlator of random 
Rashba SO interaction C(k — k'). 



diagrams in Fig.l. We neglect ladder corrections for spin 
current vertex since we assume isotropic scattering by 
impurities. On the other hand, corrections to the ver- 
tices corresponding to the random spin-orbit coupling are 
small by the parameter Rji <C 1, where R is a character- 
istic length of the fluctuations in the Rashba interaction, 
and I is the electron mean free path. This situation is 
completely different from the case of constant SO inter- 
action. Since the spin current operator does not include 
any anomalous term, there are no diagrams with the cor- 
responding vertices including the random Rashba field. 
In other words, the spin Hall effect is due to a spin- 
dependent correction to the distribution function only, 
as will be explicitly verified later in the text by consider- 
ing the kinetic equation for the density matri x 23 ' 30 . 

We assume that the field is oriented along the y-axis 
(A x = 0, A y 7^ 0) and calculate the spin-Hall conduc- 
tivity fi S H defined as = <7 s nE y . One can easily verify 
that all other components of the spin current are equal to 
zero. Calculating the diagrams, taking the trace in the 
spin subspace, and integrating over the electron energy 
s, we obtain in the static limit oj — > 0, 



CsH 



27T7TJ 



(10) 



k x (k x + k' x ) C(q) G£ (G£_ q - G£_ q ) G 



k ■ 



kq 



where C(q) = C(q) for an isotropic system. We consider 
the experimentally relevant case of weak SO coupling, 
where To <C tso and, therefore, t is very close to To. 
For the states close to the Fermi surface, the difference 



2TmG k ? _ q can be presented in the form 



r<R 

U k-q 



"k-q 



-27ri<5 (e F - £k- q ) ■ 



(11) 
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which reflects the energy conservation. By using the re- 
sulting identity, 



k 2 - (k - q) z ) /2m 
2m 



(12) 



2_ q 2 



[S (6 - Ox) + 5 (0 - 2 )} 



with denoting the angle between k and q, and 0\^ 
being two solutions of cos^ 2 = q/2k, we arrive upon 
integrating over k at 
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2 A- r 



C(q) 



q 2 dq, 



(13) 



where hp is the Fermi momentum. Taking into account 
the formula for the spin relaxation time r s due to random 
SO coupling, derived in Reb 3 -^ for £ 3> R: 



1 m 



2 A/ 



C(g) ^Ak 2 F -q 2 dq, (14) 



we obtain 



C S H = ~. • 

47T T, 



(15) 



Equation (|15|) is our main result. It shows that ct s h is 
non-universal and depends on both the disorder due to 
impurities and random SO coupling. However, it is not 
zero under very general assumptions of our model, which 
is qualitatively different from <t s h = for the uniform 
Rashba coupling^. Of course, our finite ct s h is not in 
contradiction to the result of Ref^ since, in contrast to 
Reb2£, the SO coupling is disordered here. 

The above result can also be obtained with the kinetic 
equation for random Rashba SO interactio n 23 ! 31 , which 
allows a better insight into the problem. The kinetic 
equation for the density matrix pk includes the usual 
field-dependent term eE-dpok/dk, which is responsible 
for the conductivity. Here the unperturbed density ma- 
trix is pok = Ifo(sk)i where /o(e) is the Fcrmi-Dirac dis- 
tribution function. Another source of the perturbation in 
the electron distribution under external field E is due to 
the spin-dependent scattering associated with the fluctu- 
ating Rashba SO interaction, as given by the third term 
in the right- hand-side of Eq.(J7|). As discussed above, the 
electric field produces a random field acting on electron 
spin, which can be treated in the collision integral. 

Using the matrix elements in Eqs. ([5]), ©, and assum- 
ing that the perturbation due to SO interaction is small, 
we find the following expression for the collision integral: 

St p k = - (pou - pk) + St^Jpk, (16) 

T 

St^p k = 2n £ (iW^Fk + ^ 

k' 

(pov - Pok)<5(£k - £k' + u), (17) 

where St'^pk is the contribution from the random 
Rashba field. Since we consider the linear response to 



E y , in the last term we take the equilibrium density ma- 
trix, with pok' — Pok = wdfo(s)/de. Upon the averaging 
over the SO disorder, we obtain, 

%^^ + W^V^=~G(q)(k x + k' x )a z E y . (18) 

Lu 

This term is the driving force for the spin current, as 
shown in Fig. 2. The corresponding contribution to the 
collision integral can be obtained with Eq. (IT2l as 



St 



- c d/o 1 1 ( 
de k t s \ k 



(19) 



We present the density matrix as pk = Pok + ^Pk + S]g_a z 
and with Eqs. ([TB"]) . (JTT|) . find for the steady state: 



■5k 



dfo k x t 
" 6 de Ft, y ' 



(20) 



describing spin split of the Fermi surface, corresponding 
to Fig. 2. In the stationary state the spin-dependent force 
due to the external field and random spin-orbit coupling 
in Eq. (jT5)) is balanced by the friction force due to the 
disorder, proportional to 1/r. 

Having found the distribution function, we can calcu- 
late the spin current. First, we calculate in the spin space 
Trpkia! with j* = k x a z /2m and obtain: 



rr, - z df t _k% „ 
Trp kJx = -e— v -^E y . 



The total spin current is then equal 

d 2 k 



■I 



Tr p k f x 



(2^ 



e t 

47T T, 



(21) 



(22) 



This result leads to <t s h equivalent to Eq . (|T5l) . 

To consider an example, we assume the following form 
of the correlator C(q)££±: 



C(q) = 2ir(\ 2 )R 2 



(23) 



achieved by doping quantum wells with charged impuri- 
ties. In the semiclassical limit of long-range correlations, 
fcpi? S> 1, the integral in Eq. (fT4|) becomes 



1 

t. 



2k f 



C(q)dq = 4m(\ 2 )k F R, (24) 



and the resulting spin Hall conductivity is 
a S H = — rriT (A 2 ) k F R- 



(25) 



When k F R <S 1, the system is always in the dirty limit 
of very long SO coupling-determined relaxation times. 
As a result, the spin Hall conductivity is suppressed and 
tends to zero as (kpR) 2 with decreasing kpR. This is a 
general feature of the finite-range correlators (similar to 
that in Eq. (|2"3"f ) , the effect of which vanishes due to fast 
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FIG. 2: Preferable spin-dependent scattering direction by the 
effective potential in left-hand-side of Eq. (|f 8[) . shown by hor- 
izontal arrows attached to the circles. White and filled circles 
correspond to spin components cr 2 = f and a z = —1, respec- 
tively. As a result, the Fermi line becomes spin split with the 
preferable concentration of spin-up electrons at k x > and 
spin-down ones at k x < 0, leading to the SHE. 



oscillations of the Rashba parameter on the spatial scale 
of the electron wavelength. 

Now wc can qualitatively discuss the clean limit r ^ 
Tgo- Here the spin conductivity is finite and does not de- 
pend on the magnitude of fluctuating Rashba field since 
both the relaxation rate and gain due to the external 
field are proportional to (A 2 ). However, the clean limit 
requires a separate analysis of the relaxation timescales, 



which will be considered elsewhere. 

In conclusion, we have shown that the random Rashba 
spin-orbit interaction can generate spin Hall effect, even 
in the presence of impurities. This behavior is distinct 
from that found for spatially uniform Rashba interaction, 
where in the limit of small impurity concentration, the 
potential scattering by impurities totally suppresses the 
spin Hall effect. In contrast to Ref. [33, where it was 
found that for the linear Rashba coupling this suppres- 
sion is a result of a sum rule for spin conductivity, here 
the corresponding rule cannot be established and the re- 
sulting spin current is not suppressed. We mention that 
the comparison of conventional and torque-related 2 ^ def- 
initions of spin current shows that the result in Eq.(|15p 
is independent of definition. 

In systems with nonzero spin polarization, arising for 
instance due to finite magnetization, the above discussed 
SHE is closely related to the anomalous Hall effect. If 
the concentrations of spin-up and spin-down electrons are 
different, spin separation leads to electric current which 
gives rise to the anomalous Hall effect. Since the disorder 
in spin-orbit coupling leads to the spin current, it can 
cause the anomalous Hall effect, too. 
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